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ABSTRACT 


In  a previous  publication  from  this  Iciboratory,  the  Rouse-Bueche- 
Zimm  molecular  theory  of  viscoelasticity  has  been  extended  by  using  a 
transient  network  model  to  apply  to  monodisperse  polymers  with  chain 
entanglements.  Effects  of  the  entanglements  were  modeled  both  by  the 
enhanced  frictional  coefficients  and  by  the  additional  elastic  couplings 
resulting  from  the  transient  entanglement  network.  For  binary  blends 
consisting  of  two  monodisperse  polymers  with  different  molecular 
weights,  additional  modifications  are  now  required  to  predict  their 
linear  viscoelastic  behavior.  It  is  recognized  that  entanglements 
not  only  may  form  between  chains  of  the  same  lengths,  but  also  between 
those  of  different  lengths.  For  the  latter  case,  the  longer  chain 
will  in  fact  have  the  frictional  coefficient  of  the  shorter 
chain  at  the  point  of  entanglements . The  frequency  of  such 
interactions  is  assumed  to  be  proportional  to  the  weight  ratio  of 
the  respective  component  chains  in  the  blend.  Equations  of  motion  are 
formulated  for  each  component  and  solved  numerically  for  the  relaxation 
spectra.  Linear  viscoelastic  parameters  such  as  the  dynamic 
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mechanical  moduli,  stress  relaxation  moduli,  and  zero  shear 
viscosity  can  then  be  computed  for  the  blends  by  linear 
summation  of  those  of  the  components.  The  reduced  steady 
state  shear  compliance  of  the  blends  can  also  be  computed 
from  knowledge  of  the  component  relaxation  times.  Results 
are  found  to  be  in  excellent  agreement  with  literature  data  on 
polystyrene,  poly (dimethyl  siloxane)  and  poly (methyl  methacrylate). 

INTRODUCTION 

The  molecular  theory  of  viscoelasticity  for  isolated 
polymers  expounded  by  Rouse  (1),  Bueche  (2),  and  Zimm  (3)  has 
previously  been  extended  by  a number  of  workers  (4-9)  to  predict 
the  viscoelastic  properties  of  bulk  polymers  with  molecular 
weights  exceeding  their  critical  entanglement  molecular  weights, 

M^.  In  recent  series  of  publications  from  this  laboratory  (10-12), 
several  models  were  considered  which  incorporate  the  effects 
of  enhanced  friction  and  elastic  coupling  of  the  entangled 
chains.  The  most  successful,  and  intuitively  most  reasonable, 
model  envisions  the  entangled  systems  as  forming  a three  dimen- 
sional transient  network  (12) . The  predicted  linear  visco- 
elastic properties  for  monodisperse  entangled  polymers  were 
found  to  be  in  excellent  agreement  with  literature  data  for  a 
number  of  systems. 

Attempts  to  describe  the  viscoelastic  properties  of  polymer 
blends  from  those  of  their  known  monodisperse  fractions  are 
motivated  by  its  potential  engineering  applications  to  tlie  mixed 
system.  In  addition,  an  understanding  of  the  multitudinous  inter- 
actions among  the  components  may  provide  a possible  means  toward 
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predicting  the  viscoelastic  properties  of  the  polydisperse 
polymers.  A number  of  workers  invoked  some  empirical  blending 
rules  to  obtain  the  relaxation  spectra  of  mixtures  from  their 
constituents  (13-18).  However,  polymer  molecules  should  behave 
differently  in  a blend  than  in  their  pure  state  because  of  the 
interactions  with  the  coexisting  components  of  different  chain 
lengths.  The  response  and  motion  of  the  fractions  must  therefore 
be  mutually  affected  by  all  the  components  present  in  the  blend. 
The  resulting  relaxation  spectra  associated  with  any  given 
component  must  be  different  in  a blend  from  that  in  the  pure, 
unmixed  state.  Merely  mixing  the  unmodified  relaxation  spectra 
of  the  components  to  form  the  relaxation  spectra  for  the 
blend  without  taking  into  account  of  the  aforementioned 
mutual  interactions  between  the  components  would  be  inadequate 
in  describing  the  viscoelastic  behavior  of  polymer  blends. 

In  this  work,  we  shall  apply  the  transient  network  model 
(12)  to  binary  blends  by  incorporating  the  effect  of  the  co- 
existing components  on  chain  dynamics.  We  are  then  able  to 
obtain  the  modified  relaxation  times  for  each  component.  The 
linear  viscoelastic  properties  of  the  blends  are  subsequently 
obtained  by  summing  properties  of  the  corresponding  component 
calculated  by  the  modified  relaxation  times. 

THEORY 

In  a polymer  blend  made  of  monodisperse  fractions  whose 
molecular  weights  exceed  the  critical  molecular  weight,  all 
the  polymer  molecules  entangle  with  their  surrounding  neighbosc; 
to  constitute  a three  dimensional  network  structure.  The 


difference  between  this  kind  of  network  and  the  ones  formed 
by  the  monodisperse  components  alone  is  that  in  the  blend 
entanglements  may  be  formed  from  chains  of  different  lengths. 

The  chain  dynamics  in  these  systems  may  be  expected  to  differ 
significantly  from  the  monodisperse  polymers. 

Following  the  transient  network  model  for  monodisperse  poly- 
mers, the  equation  of  motion  for  an  unentangled  bead  B is  (12) 


kT  31nij; 
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<b^> 
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(1) 


where  k^0,N,  N is  the  total  number  of  beads  in  the  chain  and 
2 

<b^>  IS  the  mean  square  end-to-end  distance  of  a statistical 

segment,  and  are  respectively  the  x-component  of  the 

coordinates  and  velocity  of  the  kth  bead,  V®  is  the  velocity 

^k 

of  the  medium  surrounding  the  kth  bead,  f_  is  the  frictional 

O 

coefficient  of  the  unentangled  bead,  and  t)*  is  the  segmental 
distribution  function. 

For  an  entangled  bead  A,  the  corresponding  equation 
assumes  the  following  form: 


k. 

3 


kT  3lnii' 


* 2x3  -3-1' 


3kT 


<b 


2> 


(x.-x^ ) 
3 1 


(2) 


where  t . . 

31 


(3) 


I 


r 

5. 

and  is  the  frictional  coefficient  of  the  entangled 

j 

bead  at  x ^ . The  fourth  term  on  the  right-hand  side  of  Eq.  2 
accounts  for  the  elastic  coupling  forces  acting  between  all 
pairs  of  entangled  beads  (12).  They  are  inversely  propor- 
tional to  the  interbead  distance  as  shown  in  Eq.  3.  o is  a 
parameter  to  account  for  the  fact  that  the  entanglement  is  not 

as  effective  as  a permanent  crosslink  in  transmitting  force,  ' 

since  entangled  chains  can  slip  by  each  other.  Its  value  is  i 

between  0 and  1.  Throughout  the  calculations,  the  enhanced 

frictional  coefficient  f^  is  assumed  to  have  the  high-friction-  , 

j ^ 

inside  distribution,  i.e.,  entangled  beads  near  the  chain  middle 
have  higher  frictional  coefficients  than  the  ones  at  the  chain 

ends.  This  type  of  distribution  is  not  only  intuitively  reason-  ! 

able,  but  also  has  been  found  to  fit  the  experimental  data  very 

1 

well  for  monodisperse  polymers  (12) . 

In  a polymer  blend  made  of  monodisperse  fractions  having 
different  molecular  weights,  further  modifications  of  the  equation 
of  motion  must  be  imposed.  For  simplicity  the  following  derivation  = 

will  be  given  for  a binary  blend.  Extension  to  multiple  blends 
can  be  similarly  analyzed.  We  begin  by  assuming  that  the  blend 
is  a homogeneous  one,  and  shows  no  evidence  of  microphase  sepa- 
ration. This  assumption  may  be  considered  valid  for  the  type  of 
blends  under  consideration  , i.e.,  blends  of  different  molecular 
weight  fractions  of  the  same  polymer.  The  absence  of  heterogeneous 
domain  formation  is  one  of  the  basic  tenets  of  the  RBZ  model. 

For  a given  component  chain  in  a binary  blend,  it  can  form  en- 
tanglements with  another  chain  of  the  same  length  or  that  of 

a different  length.  The  probability  of  entanglement  formation  | 
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between  chains  belonging  to  different  components  and  those  belonging 
to  the  scune  component  in  a well-mixed  blend  is  expected  to  be  propor- 
tional to  the  abundnace  of  the  components  if  we  assume  that  the  ten- 
dencies for  chain  segments  to  entangle  are  the  same  irrespective  of 
the  chain  lengths  as  long  as  they  exceed  the  critical  molecular  weight. 
The  ratio  of  in ter -component  and  intra-component  entanglements  is 
therefore  dependent  upon  the  composition  of  the  blends,  being  propor- 
tional to  the  weight  ratio  of  the  components  according  to  our  assumption. 
Here,  inter-component  entanglement  is  defined  as  entanglement  formed 
between  two  beads  belonging  to  chains  of  different  lengths,  and  intra- 
component entanglement  corresponds  to  that  between  chains  of  the  same 
length . 

Entanglement  formation  restricts  chain  motion  of  both  chains  at  the 
entangled  point.  The  extent  of  retardation  is  the  Scune  for  both  chains. 
At  an  inter-component  entanglement,  motion  of  the  longer  chain  is 
facilitated  because  it  is  entwined  with  a shorter  chain.  This  serves 
to  reduce  the  frictional  coefficient  of  the  longer  molecule  at  the  inter- 
component entanglement  point.  In  other  words,  the  longer  molecule  now 
has  the  smaller  frictional  coefficient  previously  associated  with  the 
shorter  molecule  while  the  frictional  coefficient  of  the  shorter  molecule 
remains  unchanged.  This  is  a reasonable  assumption,  since  the  frictional 
coefficients  must  be  the  same  for  both  beads  from  either  chain  at  the 
point  of  entanglement.  Because  of  the  greater  mobility  of  the  shorter 
chain,  the  entangled  beads  on  both  chains  now  are  less  restrained.  We 
can  now  write  the  equation  of  motion  for  such  a system  in  matrix 
notation  as: 


_ ,.s  n ..,“1  31nii)  ,,-1  „ „ 

? = Yx  ■ ? “9X  B ? ?e  - 
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where  X and  X are  column  vectors  of  velocity  and  position  of  the  beads 
respectively,  V®  is  the  x-component  of  the  velocity  of  the  surrounding 
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is  the  diffusion  coefficient,  o„  = 
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is  the  inverse  of  the  diagonal  matrix 
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where  <5^  = j 


A . 
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and  fj.  now  is  the  modified  high-f riction- 
B 

inside  distribution  of  the  enhanced  frictional  coefficients. 

For  a given  component  in  a binary  blend,  a number  of  the  6j's 
are  chosen  randomly  to  be  replaced  by  those  corresponding  to 
the  other  component  so  as  to  represent  inter-component  entan- 
glements. The  frequency  of  such  exchange  is  proportional  to 
the  amount  of  the  latter  component  present  in  the  blend.  is 

the  modified  version  of  the  nearest-neighbor  matrix  for  an  en- 
tangled polymer.  As  an  example,  the  form  of  for  a chain  with 
alternating  entangled  and  unentangled  beads  is: 


(6) 
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whore 


N 

j = l 

(i-1)  , i , (i+l) 


(7) 


Wo  can  write  a matrix  equation  for  each  component  in  the 
bieni.1  and,  upon  solving  it,  obtain  the  relaxation  times  of  the 
component  qiven  by 
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where  > is  the  p-th  eiqenvalue  of  D ^ Z • The  interactions 
between  components  have  been  incorporated  into  the  matrix  equa- 
tions, the  resulting  component  relaxation  times  are  therefore 
the  modified  relaxation  times.  The  degree  of  modification 
depends  on  the  composition  .is  explained  before,  the  calculated 
relaxation  times  for  the  components  hold  true  only  for  the 
particular  composition  for  which  the  calculations  have  been 
carried  out.  The  rela.xation  times  for  the  same  component  polymer 
chain  vary  with  the  composition  of  the  blend. 

Further  refinement  must  be  introduced  to  our  basic  theore- 
tical framework  upon  recognition  of  the  fact  that  these  inter- 
component entanglements  can  take  place  at  different  positions 
along  a given  chain.  Since  we  have  adopted  the  high-friction- 
inside  distribution  as  our  basis,  replacing  frictional  coeffi- 
cients means  utilizing  a different  set  of  numbers  when  the  en- 
tanglements occur  at  different  positions  along  both  chains. 

Two  approaches  can  be  employed  to  tackle  this  problem. 
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The  first  method  is  to  select  randomly  the  beads  whose 
frictional  coefficients  are  to  be  changed,  and  then  calculate  the 
relaxation  times  from  which  the  various  viscoelastic  properties 
can  be  found.  The  whole  procedure  is  repeated  many  times  and  the 

i 

resulting  viscoelastic  properties  are  averaged  in  an  attempt  to 
produce  a good  representation  of  the  system.  This  is  the  "post- 
averaging" method. 

The  second  method  assumes  that  all  the  entangled  beads  along 
I a given  chain  have  equal  probability  of  forming  intercomponent 

' entanglements  with  any  of  the  entangled  beads  of  the  other  components. 

i 

As  an  example,  consider  a 50/50  weight  percent  binary  blend.  A 
i given  entangled  bead  on  a given  chain  is  involved  in  inter-component 

entanglement  half  of  the  time  and  intra-component  entanglement  the 
other  half  of  the  time  over  a long  period  of  observation. 

! 

i Equivalently,  for  a mixture  containing  a huge  number  of  molecules, 

for  the  same  entangled  point,  half  of  the  molecules  are  involved 
f in  inter-component  entanglement  formation  and  the  other  half  in  intra- 

component entanglement  at  any  given  instant.  We  can  therefore 
"pre-average"  the  frictional  coefficients  and  use  the  resulting 
j average  to  set  up  the  matrix  equations  as  a representation  of  the 

i collection  of  chains  in  the  mixture. 

i 

I 

* Suppose  that  our  example  of  a 50/50  (by  weight)  binary  1 

blend  is  made  of  components  having  10  and  20  entanglements 
respectively.  Based  on  the  "equal  ’entangling  tendency"  assiimption, 
five  out  of  the  ten  entanglements  in  the  shorter  component  are 
involved  in  in ter -component  entanglements  on  the  average,  and 


ton  out  ot  twenty  toi  the  lotKiet  oomiHMient.  The  post -aver uij i inj 
method  randomly  selects  ten  of  the  matrix  correspond- 

in>.j  to  the  lofujer  component  anii  lep  laces  them  by  those  coi' respond  i luj 
to  the  shorter  component.  The  uuunted  values  are  also  picked 
randomly  from  the  '■'j''’  *^1  i^l'orter  oi.miponents . The  result  inij 

viscoelastic  parameters  are  calculated  aiul  stoted  in  the  memory 
bank,  and  at  the  end  of  the  computation,  tiu'  stored  data  ate  avei.uied 
to  obtain  the  "post-avera-.jed"  propeitie^. 

It  the  random  selection  process  is  repeated  many  times,  event- 
ually any  ontanylement  aloiu)  the  longer  chain  will  assume  its  ori«.iinal 
frictional  coefficient  half  of  the  times,  while  tor  the  other  half 
of  the  times  it  will  have  frictional  coefficients  belon^.! i nv-i  to  the 

shorter  component.  Since  we  adopted  the  hitjh-f  riot  ion-inside  dis- 

I 4 

tribution  (12),  i.e.  f^  = 1 , wliere  I is  the  position  of  the 

J 

entanglement  counted  from  chain  ends,  the  proposed  preavera^j i luj 
method  assiijns  the  frictional  coefficients  accordiiuj  to  Uie  followin«.j 


rule : 


f . = W . 1 + W. 

Aj  J k 


when  N is  even;  and 


^ 4 

V IN 

e,  ■ 


f = W l"*!-  WUN^  2 + !)■’  Y.  L"*)  ] /N, 

^ ^ h=l  "k 


when  N is  odd.  In  eos . 9 and  10,  W is  the  weioht  fraction  and  N 
^k  ’ 

is  the  number  of  entant.j lements  of  the  loiujer  component  while  sub- 
script k refers  to  the  shorter  comv'onent  in  the  blend.  Fq . 4 ensures 


that  the  distribution  is  still  h iqh-f r ict ion  inside  while  the  effect 
of  inter-component  interactions  is  adequately  taken  into  account . 
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Solution  oi  the  pos t -uvei  i no  method  is  time-consunung  and 

expensive  in  the  cost  of  computa t ions . One  case  has  been  tried  and 
was  found  to  approach  the  pre-avei  ay iiu7  result  as  tJie  number  of 
reiteration  of  the  random  selection  process  increases.  The  pre- 
averayiny  method  is  thus  chosen  for  the  subsequent  calculations. 

Once  the  modified  relaxation  times  are  found,  they  can  be 
used  to  calculate  all  the  linear  viscoelastic  properties  such  as 
relaxation  modulus,  dynamic  moduli,  steady-state  shear  compliance 
and  zero-shear  viscosity  by  the  well-known  relations  (10).  The 
blend  properties,  e.\cept  for  the  steady-state  shear  compliance, 
are  iust  the  linear  sums  of  the  corresponding  properties  of  compon- 
ents because  the  relaxation  times  for  each  component  are  obtained 
from  the  already  modified  matrix  equation  to  meet  conditions 
e.xistiny  at  the  particular  composition  considered.  The  steady-state 
shear  compliance  of  the  blend  is  obtained  from  the  resulting  relax- 
ation times  of  the  components  by  the  appropriate  equation  to  be  shown 
later  . 
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RESULTS  AND  DISCUSSION 


In  comparing  the  theoretical  prediction  with  experimental 
data,  one  of  the  computed  curves  was  shifted  in  coordinates 
until  the  terminal  regions  overlapped  with  data.  All  sub- 
sequent curves  were  not  adjusted  further  once  the  fitting 
parameters  have  been  chosen.  The  dynamic  storage  and  loss 
shear  moduli  were  computed  by: 
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Figure  1 compares  the  theoretical  results  with  data  on 
polystyrene  blends.  The  solid  curves  are  the  computed  curves. 
The  symbols  represent  experimental  data  by  Onogi  et  al.  (19) . 

To  calculate  the  number  of  entanglements  for  a given  mono- 

4 

disperse  system  (N^) , is  assumed  to  be  1.7x10  (10-12). 

The  two  curves  on  the  outside  are  for  the  pure  components 

4 5 

of  molecular  weights  of  5.9x10  and  6.16x10  respectively. 

We  can  see  that  the  experimental  data  on  the  shorter  component 
have  not  been  accurately  predicted.  The  reason  is  that  the 
equivalent  number  of  entanglements  for  this  component  should 
be  about  3.5.  However,  the  computer  program  can  only  calculate 
chains  with  integral  number  of  entanglements  and  we  are  forced 
to  use  = 3.  The  middle  curve  is  a 50/50  (by  weight)  blend 
of  the  two  components.  The  two-plateau  tendency  in  the  rubbery 
region  of  the  dynamic  storage  moduli  is  successfully  predicted. 


In  addition,  the  coirospondiiuj  n\aximnn\  near  losl  »•  0 is  also 
in  oxcollent  acjroomont  with  data.  Those  phenomena  are 
frequently  observed  for  binary  bUsiuis  of  two  components  haviiuj 
very  different  molecular  weiqlits.  I'iq.  2 shows  results  on 
another  sot  of  polystyrene  binary  blends  by  the  same  authors  (19). 
The  blends  are  20/80,  40/00,  00/40,  80/20  l)y  wt'iqht.  Here  the 
shorter  component  has  a molecular  wi'iqht  cor respondinq  to  2.5 
entanqlements.  It  is  assumed  to  have  2 ent anq Icments  in  our 
calculations.  Aqreement  betwi'on  theory  aiui  experiment  is  »\qain 
very  qood. 

I'iqure  3 compares  the  calculated  curves  with  a set  of 

dynamic  mechanical  data  on  blends  of  poly (methyl  methacrylate) 

3 4 

of  molecular  weiqhts  4.1x10  anil  1.74x10  . The  blends  as  in- 
dicated by  the  middle  curves  are  90/10,  80/20,  80/40,  30/70 
percent  by  weiqht,  the  first  number  in  the  p.iir  of  numbers  is 
the  weiqht  fraction  pertaininq  to  the  liqhter  component. 
for  PMMA  is  assumed  to  be  7,000  (21).  Aqreement  between  theory 
and  experiments  is  qood  except  that  the  experimental  data  exhibit 
a more  qradual  terminal  reqion,  possibly  dvie  to  its  broader 
heteroqeneity  index  (11.1.  = 1.25). 

The  stress  relaxation  moduli  calculated  by  the  well  known 
relation 


R 


(14) 


are  shown  in  Fiqure  4.  Also  qiven  are  the  experimental  data 
of  Murakami  and  Ono  (21)  tot  monodisperse  polystyrene  fractions 


b 5 

with  molecular  weights  ot  1.87x10  aiul  5.8Jxl0  , as  well  as  their 
72/25  and  50/50  blends.  The  agreement  between  theory  and  experiment 
is  good,  except  that  the  small  second  plateaus  were  not  appartuit 
in  the  predicted  curves. 

To  compute  the  steady-state  shear  compliance,  we  lecall  that 
by  definition 

= NkT);  1 "/ (NkT.\  1 . ) (15) 

e r i 

tor  a monodisperse  polymer.  For  a binary  blend,  then 
N '■  1 ^ + N '■  1 ^ 

jO  ^ ^ (H,) 

RT(N^V  1 . + N^y.  i j)  *■ 

where  i^'s  and  relaxation  times  of  components  1 and  2 

respectively,  and  Nj,N2  are  numbers  of  molecules  of  components  1 
and  2 per  unit  volume.  In  other  winds, 

jO  = (17) 

^ NRT(f  i . + f ^y  i . ) ‘■ 

where  f^^  and  f^  are  mole  fractions  of  components  1 and  2,  i.e., 

W,/M 

f,  = - — — T—  and  f.,  = 1-f,.  The  reduced  steady-state  shear 

1 Wj^/Mj^+W^/M^  2 1 

compliance  for  a binary  blend  is  thus 

t 'll?  + f.,yi'^ 

. 112]  /io\ 
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Figure  5 shows  the  Log  J - W_  data  for  three  sets  of  binary 
blends  of  polystyrene.  In  addition  to  the  data  of  Onogi  (19)  and 
Akovali  (22) , the  more  recent  data  of  Porter  and  Prest  (24)  are  also 
included.  Both  the  data  of  Onogi  and  those  of  Porter  and  Prest  show 
pronounced  maxima,  and  are  in  fairly  good  agreement  with  our  com- 
puted curves.  The  molecular  weights  in  the  binary  blend  of  Onogi 
are  4.69x10^,  1.67x10^  and  those  of  Porter/Prest  are  9.72x10^, 

4.]  1x10^.  The  molecular  weights  in  the  A)covali  blend,  on  the  other  hand, 
5 5 

are  1.25x10  , 2.67x10  . The  molecular  weight  ratio  is  much  lower 
than  those  in  the  former  two  blends.  The  J now  exhibit  a much 
gentler  maximum  when  plotted  against  The  computed  curves  are 

still  in  reasonable  agreement  with  the  data,  although  the  predicted 
maximum  is  sharper.  Of  course,  if  the  molecular  weight  ratios  in 
the  binary  blends  decreases  further  toward  unity,  then  a horizontal 
straight  line  would  eventually  be  reached  (a  monodisperse  system) . 

Thus  the  trend  is  consistent  with  expectation.  In  Figure  6 we  com- 
pare similar  data  of  poly (methyl  methacrylate),  which  are  seen  to 
be  in  excellent  agreement  with  the  theoretical  curve. 

The  zero-shear  viscosity  can  be  computed  from  the  relaxation 


spectra  by 


^o  ~ N 

e p 


The  computed  values  are  then  plotted  in  Figure  7 against  the  weight 
average  molecular  weight  of  the  blends  containing  different  amounts  of 
the  two  component  monodisperse  polystyrene  fractions.  The  predicted 
curve  is  a straight  line  with  a slope  of  3.70,  which  is  in  good 
agreement  with  the  experimental  data  of  polystyrene  represented  by 


the  circles  (20)  . 


T 
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The  variation  of  the  blend  viscosity  as  a function  of  com- 
position is  an  important  engineer! ny  property  and  also  a good 
criterion  for  testing  our  theory.  Friedman  and  Porter  (23)  reviewed 
an  exhaustive  list  of  mixing  rules,  and  compared  them  with  experi- 
mental data.  Following  their  convention,  we  plot  in  Figure  8 the 
viscosities  of  binary  polystyrene  blends  against  the  weight  fraction 
of  the  heavier  component.  The  circles  are  experimental  data  (19) 
and  the  solid  curves  are  the  theoretical  predictions.  The  numbers 
in  the  parenthesis  give  tl\e  equivalent  entanglement  bead  number 
for  the  two  components  from  which  the  blends  are  made  of . 

Agreement  between  theory  and  experiments  is  seen  to  be  excellent. 
Figures  9 and  10  show  similar  plots  for  two  other  polymer  systems, 
i.e.,  poly(methyl  methacrylate) ( 21 ) and  poly (dimethyl  siloxane)  (26). 
Agreement  is  again  very  good . 

For  the  cases  examined  above,  all  the  viscosity  - composition  c 
curves  concave  downwards.  However,  the  polystyrene  blends  studied  by 
A)covali  (22)  revealed  a somewhat  sigmoidally-shaped  curve.  The  only 
significant  difference  between  this  blend  and  the  others  is  that 
the  two  components  have  a much  smaller  chain  length  difference  in 
the  former  (Figure  11)  than  in  the  latter  (Figures  8-10) . Interest- 
ingly this  trend  is  predicted  by  our  theory  (Figure  11) . To  further 
scrutinize  the  effect  of  chain  length  differences  between  components 
on  the  shape  of  the  viscosity-composition  curves,  we  show  in  Figure  12 
the  computed  curves  for  a series  of  binary  blends  with  varying 
chain-length  differences.  We  note  that  if  the  lengths  of  the  two 
components  are  not  significantly  different  (up  to  a factor  of  two), 
the  curves  are  sigmoidal.  The  theory  in  fact  predicts  a minimum 
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in  viscosity  if  the  chains  are  of  comparable  lengths,  which  has 
not  been  experimentally  verified.  The  concave  downward  appear- 
ance, on  the  other  hand,  will  be  observed  for  components  whose 
chain  lengths  differ  more  than  a factor  of  two. 

In  conclusion,  we  can  state  that  the  transient  network  model 
is  successful  in  predicting  the  linear  viscoelastic  behavior  of 
monodisperse  polymers  (12)  as  well  as  their  binary  blends.  The 
important  feature  of  the  model  for  blends  is  that  it  is  not 
sufficient  to  just  use  some  kind  of  mixing  rules  for  the  pure 
components  alone.  Simple  intuitive  reflection  of  the  binary  system 
will  indicate  that  the  cross  entanglements  between  chains  of 
different  lengths  must  also  be  taken  into  account.  The  good  agree- 
ment between  the  theoretical  predictions  and  the  literature  data 
appears  to  vindicate  the  validity  of  this  model. 
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Table  1:  Parameters  Used  in  the  Calculations  and  the 

Resulting  J „ and  n • 
eR  o 


Polymer 

N+1 

Ne 

6^ 

3 

a 

“^eR 

"o 

Polystyrene 

7 

3 

0.2 

0.333 

0.972  ^ 

50/50 

blend 

II 

II 

- 

1.92x10 

69 

34 

II 

II 

0.101 

7.69x10^ 

Polystyrene 

7 

2 

0.2 

0.409 

0.508 

80/20 

blend 

II 

If 

- 

3.52 

60/40 

II 

II 

II 

- 

1.13x10 

40/60 

II 

M 

II 

- 

2.37x10 

20/50 

II 

II 

II 

- 

4.09x10 

21 

lo’ 

II 

II 

0.269 

6.29x10 

Polystyrene 

23 

11 

0.2 

8.50xl0_ 

75/25 

blend 

II 

II 

- 

4.71x10 

50/50 

II 

II 

II 

- 

1.70x10^ 

69 

34 

II 

II 

— 

6.71x10^ 

Polystyrene 

7 

2 

0.2 

0.402 

0.499 

90/10 

blend 

M 

II 

3.759 

1.09 

80/20 

It 

If 

II 

3.824 

2.82 

70/30 

II 

II 

II 

2.842 

5.66 

60/40 

II 

II 

II 

2.073 

9.62 

40/60 

II 

II 

II 

1.115 

2.10x10 

20/80 

tl 

II 

II 

0.591 

3.68x10 

21 

10 

II 

M 

0.265 

5.72x10 

Polystyrene 

15 

7 

I'* 

0.2 

0.336 

1.52x10 

90/10 

blend 

If 

II 

0.828 

1.74x10 

80/20 

II 

II 

II 

1.023 

2.51x10 

70/30 

II 

II 

II 

0.871 

3.83x10 

60/40 

II 

II 

II 

0.692 

5.67x10 

40/60 

II 

II 

II 

0.445 

l.lOxlO^ 

20/80 

II 

If 

II 

0.291 

1.85x102 

31 

15 

II 

II 

0.193 

2.82x102 

Polystyrene 

13 

6 

0.2 

0.396 

8.4 

90/10 

blend 

II 

II 

2.313 

2.58x10 

80/20 

H 

II 

II 

1.808 

7.76x10 

70/30 

H 

II 

It 

1.246 

1.64x102 

60/40 

II 

II 

II 

0.874 

2.85x10^ 

40/60 

N 

N 

•1 

0.474 

6.27x10!; 

20/80 

M 

II 

•1 

0.261 

1.11x10^ 

49 

24 

tl 

If 

0.131 

1.73x10^ 
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Table  1 (cont'd) 
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Poly (methyl- 
Methacrylate) 


13  6 

90/10  blend 
80/20 
70/30 
60/40 
40/80 
20/80 

51  25 


tl 


II 


0.1  0.410  1.04x10 

" 1.369  4.03x10 

" 2.143  1.29x10 

" 1.428  2.76x10 

" 0.994  4.82x10 

" 0.532  1.07x10 

" 0.289  1.90x10 

••  0.127  2.05x10 


Poly (dimethyl- 

7 

3 

I" 

Siloxane) 

80/20 

blend 

II 

60/40 

II 

II 

40/20 

If 

II 

20/80 

II 

II 

79 

39 

II 

0.972 

1.02x10 

6 . 61x10 

2.44x10 

6.72x10 

9.75x10 
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Figure  1 


Figure  2 


Figure  3 
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CAPTIONS  OF  THE  FIGURES 

Plots  of  dynamic  storage  moduli  G'  (i»J)  and  loss 
moduli  G"  (ii>)  vs  frequency  for  polystyrene.  The 
solid  curves  are  computed  curves.  The  symbols  are 
experimental  data  (19) . The  two  outer  curves  are 
for  the  pure  components.  The  numerals  at  the  bottom 
of  the  curves  indicate  the  numbers  of  entangled  beads. 

The  curve  in  the  middle  (triangles)  is  a 50/50  (by 
weight)  blend. 

Plots  of  dynamic  storage  moduli  G' (w)  and  loss 
moduli  G" (w)  vs  frequency  for  polystyrene.  The 
solid  curves  are  computed  curves.  The  symbols  are 
experimental  data  (19).  The  two  outer  curves  are 
for  the  pure  components.  The  numerals  at  the  bottom 
of  the  curves  indicate  the  numbers  of  entangled  beads. 

The  curves  between  those  for  pure  components  are 
for  20/80,  40/60,  60/40,  80/20  (by  weight)  blends. 

Plots  of  dynamic  storage  moduli  G' (w)  and  loss 
moduli  G" (w)  vs  frequency  for  poly (methyl  methacrylate). 
The  solid  curves  are  computed  curves.  The  symbols  are 
the  experimental  data  (20) . The  outer  two  curves  are 
the  pure  components.  The  numerals  at  the  bottom  of 
the  curves  indicate  the  number s of  entangled  beads. 

The  curves  between  those  of  the  pure  components  are 
for  90/10,  80/20,  60/40,  30/70  (by  weight)  blends. 
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Figure  4.  Stress  relaxation  moduli  for  polystyrene.  The 
solid  curves  are  computed  curves.  The  symbols 
are  experimental  data  (22)  . The  two  outer  curves 
are  for  the  pure  components.  The  ones  between  the 
pure  component  curves  are  for  75/25,  50/50  (by 
weight)  blends. 

Figure  5.  Log-log  plot  of  the  reduced  steady-state  shear  com- 
pliance of  polystyrene  vs  weight  fraction  of  the 
heavier  component.  Triangles,  circles  and  squares  are 
experimental  data  by  Porter  (26),  Onogi  (19)  and 

Akovali  (22)  respectively.  Curves  were  all  computed 
by  theory . 

Figure  6.  Log-log  plots  of  the  reduced  steady-state  shear 

compliance  of  poly (methyl  methacrylate)  vs.  weight 
fraction  of  the  heavier  component.  Circles  are 
data  of  Onogi  (20) . The  curve  was  calculated  by 
theory. 


Figure  7.  Log-log  plots  of  zero-shear-rate  viscosity  vs 
molecular  weight.  The  lines  were  the  computed 
curves.  The  circles  are  the  experimental  data  (19) 
for  pure  polystyrene  components  and  blends,  all  of 
which  fall  on  the  same  straight  line  with  a slope  of 
3.7.  The  triangles  are  data  on  poly (methyl  metha- 
crylates) (20)  . 


Figure  8.  Zero-shear  viscosities  of  two  polystyrene  binary  blend 


systems  as  a function  of  the  weight  fraction  of  the 
heavier  component.  Circles  are  experimental  data  of 
Onogi  (19) . Solid  curves  are  theoretical  predictions. 
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Figure  9.  Zero-shear  viscosities  of  poly (methyl  methacrylate) 
binary  blends  as  a function  of  the  weight  fraction 
of  the  heavier  component.  Circles  are  experimental 
data  of  Onogi  (21) . Solid  curves  are  theoretical 
predictions. 

1 

Figure  10.  Zero-shear  viscosities  of  poly (dimethyl  siloxane) 
binary  blends  as  a function  of  the  weight  fraction 
of  the  heavier  component.  Circles  are  experimental 

data  of  Prest  (25) . Solid  curves  are  theoretical  j 

I 

predictions.  j 

i 

Figure  11.  Zero-shear  viscosities  of  polystyrene  binary  blends  j 

as  a function  of  the  weight  fraction  of  the  heavier  ■ 

comp\,nent.  Circles  are  experimental  data  of  Al^ovali 
(22).  Solid  curves  are  theoretical  predictions. 


Figure  12.  Zero-shear  viscosities  of  binary  blends  for  five 
different  systems  with  different  molecular  weight 
ratios  between  the  two  constituent  components.  The 
curves  are  all  calculated  from  the  theory. 
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